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We develop a numerical tool for the fast computation of the temperature and polarization power 
spectra generated by domain wall networks, by extending the publicly available CMBACT code — 
that calculates the CMB signatures generated by active sources — to also describe domain wall 
networks. In order to achieve this, we adapt the Unconnected Segment model for cosmic strings 
to also describe domain wall networks, and use it to model the energy-momentum of domain wall 
networks throughout their cosmological history. We use this new tool to compute and study the 
TT, EE, TE and BB power spectra generated by standard domain wall networks, and derive a 
conservative constraint on the energy scale of the domain wall-forming phase transition of t] < 
0.92 MeV (which is a slight improvement over the original Zel’dovich bound of 1 MeV). 


I. INTRODUCTION 

Domain walls are two-dimensional topological defects 
that are formed when a discrete symmetry is sponta¬ 
neously broken in a phase transition [l| . The production 
of domain wall networks as remnants of phase transi¬ 
tions in the early universe is predicted in several grand 
unified scenarios. They are, however, often overlooked 
in cosmology since the energy scale of the phase transi¬ 
tion that originates the domain walls is restricted to be 
smaller than 1 MeV 0. Nonetheless, there is still room 
for domain walls to play a relevant role in cosmology: for 
instance, they have been sugrested as a possible signifi¬ 
cant dark energy contributor [3|— albeit they have since 
been ruled out as a major dark energy component [3-@| 
- and as a possible explanation for the spatial variations 
of the fine structure constant hinted at by HIRES/Keck 
and VLT/UVES data 0-0. 

Standard domain wall networks are expected to persist 
throughout cosmological history, and their presence at 
late times would necessarily leave imprints on the Cosmic 
Microwave Background (CMB). Although current CMB 
observations seem to be consistent with the inflationary 
paradigm, in which the fluctuations are seeded in the very 
early universe, they also allow for a subdominant topo¬ 
logical defect contribution 0. Domain walls, as cosmic 
strings, source metric perturbations actively throughout 
the cosmological history. For this reason, their CMB sig¬ 
natures are expected to be fundamentally different from 
those due to primordial fluctuations. In particular, cos¬ 
mic defects are expected to generate a significant vector 
component that would not be present in inflation-seeded 
scenarios, because vector modes decay rapidly in the ab¬ 
sence of a source [ll|. The B-mode polarization signal 
originated by topological defects has then contributions 


from both tensor and vector modes, and may, for this 
reason, produce an observationally relevant signal in this 
channel, despite providing only subdominant contribu¬ 
tions to the temperature and E-mode power spectra. The 
B-mode polarization channel offers thus a relevant obser¬ 
vational window for the detection of topological defects. 

Although some preliminary studies of the CMB signa¬ 
tures generated by domain walls have been done in the 
past iEMl), detailed studies of the anisotropies gen¬ 
erated by these networks are yet to be performed. In this 
paper, we develop a numerical tool to compute the angu¬ 
lar power spectrum generated by domain wall networks 
(in both the temperature and polarization channels), by 
adapting and extending the publicly available CMBACT 
code |15l - [l^ in order to accommodate these networks. 

This paper is organized as follows. In Sec. im we briefly 
review the dynamics of infinitely thin and featureless do¬ 
main walls. In Sec. Ell we adapt the Unconnected Seg¬ 
ment Model (USM) to describe domain wall networks. 
We compute the energy-momentum tensor of a domain 
wall section in Sec. IlYl In Sec. El we present and discuss 
the temperature and polarization power spectra obtained 
using this method. We then conclude in Sec. I VII 


II. DOMAIN WALL DYNAMICS 

The equation of motion of a domain wall is determined 
by the underlying field theory. However, if its thickness is 
negligible when compared to its curvature radii and if the 
domain wall is featureless — as is often the case in cosmo¬ 
logical scenarios —, the domain wall may be treated as 
a 2-dimensional surface. While moving in spacetime, an 
infinitely thin and featureless domain wall sweeps an ef¬ 
fectively 3-dimensional worldvolume. The world history 
of a domain wall may thus be represented by 
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where is a timelike parameter and and ^2 are space¬ 
like parameters. These parameters may be regarded, at 
least locally, as coordinates on the worldvolume. The 
action of such a domain wall is given by [ij 
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where a is the mass per unit area of the domain wall, 
h = det(/iab), and 


e = S'Hex^ , 


where 
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is the metric induced on the worldvolume by the back¬ 
ground (or pull-back metric). 

By varying the Nambu-Goto action in Eq. with re¬ 
spect to , one obtains the following equation of motion 
for the dynamical variables 


is the coordinate energy per unit area, and T-L = a/a. 

III. MODELING THE DOMAIN WALL 
NETWORK 


^ ^ = 0. (4) 

Let us now consider a domain wall in a flat Friedmann- 
Robertson-Walker (FRW) universe, whose line element is 
given by 

ds^ = a^{rf) {—drf' + dr"^ + -|- sin^ Od'^cj))] i (5) 

where dt] = dt/a(t) is the conformal time, t is the physical 
time, a is the cosmological scale factor, and (r, d, (f>) are 
spherical coordinates. 

The domain wall equation of motion should be invari¬ 
ant under re-parameterization of the worldvolume. In 
an expanding background, the temporal-tranverse gauge, 
with 

a;° = 77 , and x • x = 0 , (6) 

would be a natural choice. Note that, since domain walls 
are assumed to be featureless, their physical velocity is 
purely orthogonal. Let us choose, for simplicity, a set of 
spacial parameters for the domain wall that are also or¬ 
thogonal in the vicinity of the point under consideration, 
so that 

x'(i).x'W=0. (7) 

Here, x^ = (77, x), dots and '(i) represent derivatives with 
respect to conformal time and the 7-th spatial parameter 
of the worldvolume, respectively, and the latin indices in 
parentheses refer to the spatial parameters of the wall. 
In this gauge, Eq. O yields 


The main objective of this article is to devise a numer¬ 
ical tool to compute the temperature and polarization 
power spectra generated by domain wall networks. In 
order to attain this goal, different approaches might be 
followed (as was done for cosmic strings): one may ei¬ 
ther use numerical simulations of domain wall networks, 
or use a phenomenological model to describe the net¬ 
work dynamics. Here, we choose to follow the latter 
option, given the adaptability of this model based ap¬ 
proach: it allows for the description of different scenarios 
(by calibration of its free parameters), and it is not as 
constrained in terms of dynamical range as simulations 
often are. 

For cosmic strings, this approach lead to the develop¬ 
ment of a publicly available tool for the computation of 
the temperature and polarization power spectrum gen¬ 
erated by these networks — the CMBACT code [1^, [I 3 
— that is based on the phenomenological Unconnected 
Segment Model (USM) [isl - l^ to describe the energy- 
momentum tensor of the cosmic string networks. This 
tool was proven successful in predicting the shape of the 
angular power spectrum generated by cosmic string net¬ 
works obtained using other models (see e. g. [2ll,[^). It 
not only provides a good fit to the angular power spec¬ 
trum predicted using Nambu-Goto simulations, but may 
also be calibrated in such a way as to mimic that ob¬ 
tained in Abelian-Higgs simulations (by calibrating it to 
match the field theory unequal time correlators predicted 
by field theory simulations) — the so-called Abelian- 
Higgs mimic models described in [l^. Moreover, the 
bispectrum of cosmic string induced matter fluctuations 
predicted by CMBAGT is in good a gree ment with that 
computed using perturbation theory [^. The USM and 
GMBACT are, then, a flexible and robust framework to 
compute GMB anisotropies generated by cosmic strings, 
and, for this reason, we choose to extend them for domain 
wall networks. 
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A. The Velocity-Dependent One-Scale Model for 
Domain Walls 


Let us now consider a domain wall network in a FRW 
universe, and assume that the network is statistically ho¬ 
mogeneous on sufficiently large scales. In this case, two 
variables are sufficient to describe the large-scale cosmo¬ 
logical evolution of the domain wall network. One such 
variable is the root-mean-squared (RMS) velocity of the 
network, v = \J{v"^), defined as 


-2 _ Jv^ed'^^ 


( 11 ) 


where v = |x| is the microscopic velocity. The other 
dynamical variable is the characteristic lengthscale of the 
network, L, defined as 


P=^, ( 12 ) 

where p is the average domain wall energy density. 

The Velocity-dependent One-Scale (VOS) model pro¬ 
vides a quantitative description of the large-scale dynam¬ 
ics of domain wall networks, by following the cosmologi¬ 
cal evolution of these two variables. The evolution equa¬ 
tion for the RMS velocity may be obtained by averaging 
the microscopic equation of motion for a domain wall in 
Eq. ([5]). On the other hand, one may obtain an equa¬ 
tion of motion for the characteristic lengthscale of the 
network by differentiating the total energy density in do¬ 
main walls — given by 


E = aa{r]) 




(13) 


— and using Eqs. ®-(IIlD. These equations assume the 
form: 


= (1 - 
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and were derived in [23 - 1^ . (See also Refs, [l^, for 
a unified framework for the description of topological de¬ 
fects of arbitrary dimensionality.) Here, and for the rest 
of this paper, we choose to work with the comoving char¬ 
acteristic length, I = L/a. Moreover, we have introduced 
the phenomenological parameters c and n that quantify, 
respectively, the energy loss caused by domain wall inter¬ 
actions and the effect of wall curvature on their dynamics. 
In this paper, we will assume these parameters take the 
values 


c = 0.34 and /c = 0.98, (16) 

as indicated by the latest calibration of the VOS model 
against field theory simulations of standard domain wall 


networks (^ . during both the radiation and matter eras. 
Note that non-standard domain wall networks, such as 
networks with junctions, may also be describe by this 
model [^. In that case, however, these parameters would 
need to be recalibrated. 


B. An Unconnected Section Model for Domain 
Walls 

The VOS model merely provides a description of the 
large scale dynamics of a domain wall network. However, 
this is not sufficient to compute the cosmic microwave 
signature generated by domain walls: one also needs to 
characterize the energy-momentum tensor of the network 
in order to do so. To achieve this, we follow closely the 
USM for cosmic strings (Tsl - I^ . and adapt it in order 
to describe domain walls. We preserve the essential el¬ 
ements of this framework and, therefore, we will only 
briefly review its essential features. In this model, the 
network of domain walls is represented by a collection 
of uncorrelated, flat and square domain walls, that have 
been produced simultaneously at some early time. The 
VOS model is used to set the comoving length of the 
network at any given time — so that each wall has a 
comoving area P{t) —, and to fix the magnitude of the 
velocity of the segments. 

The positions of the domain walls are drawn from an 
uniform distribution in space. Moreover, the direction of 
their velocity is chosen from a uniform distribution on a 
two sphere and, since we are assuming that the domain 
walls are featureless and, consequently, that their veloci¬ 
ties are purely orthogonal, this direction also determines 
the orientation of the domain wall. 

Throughout cosmic history, a fraction of domain walls 
decays at each epoch, so that the energy loss caused by 
domain wall interactions is taken into account. As is the 
case with the USM for cosmic strings, all domain walls 
that decay at the same time are consolidated in a sin¬ 
gle section. The number of sections that decay between 
instants pi-i and rji is given by 


M{r]i)=V[n{r]^-i)-n{r]i)\ , (17) 

where V is the simulation volume, and n ( 77 ) is the number 
density of domain walls 


( 18 ) 

Here C (rj) is determined by requiring that the number of 
walls is given by V /l^ (jf) at any given time. In this way, 
at any given epoch, the number density of domain wall 
sections in the USM is in agreement with that predicted 
by the VOS model. 

The decay of domain walls must necessarily be accom¬ 
panied by a turning-off of the energy-momentum of that 







fraction. As is the case for the USM for cosmic strings, 
this effect is included by expressing the Fourier transform 
of the total energy-momentum as 
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N 


0^.(k,ry)=^[Ar(,7,)]5 0;,T°ff(r;,ry„Lf) , (19) 

i=l 

where the sum runs over all the N domain wall sections, 
0(jj^ is the Fourier transform of the energy-momentum 
tensor of the i-th domain wall section. Furthermore, 
T°® (ry, r]i, Lf) is a function that controls domain wall de¬ 
cay and turns off their energy-momentum contribution; 


r°®(ry,ry,,Lf) = 


where 


_ r, HLfVi)/v _ , 
^ In(Lf) 


( 21 ) 


and Lf is a parameter that controls the speed of domain 
wall decay. 


1 , for ry < Lfrji 

5 + 3 “ ^y) ’ <v<vi, 

0 , ioT rji < r] 

( 20 ) 


IV. THE ENERGY-MOMENTUM TENSOR OF 
A DOMAIN WALL 


With the USM for domain wall networks set up, the 
only other ingredient missing is the energy-momentum 
tensor of the domain wall sections. The energy- 
momentum tensor of an infinitely thin and featureless 
domain wall may be obtained by varying the action in 
Eq. ([2]) with respect to the metric tensor 


( 22 ) 

In the temporal-transverse gauge, it assumes the form 


(ry, x) = cr J dr](f [x'^ — {f], , ^2)] i''— ^ ^ a; -I-x (23) 


where (a;) is the 4-dimensional Dirac delta function. transform of the energy-momentum tensor of a domain 
Returning to the USM for domain walls, the Fourier wall section is, then, given by 


0^i,(k,7y) = J (ry,^i,^2) 


All 


/ : 


1-112 J 

- 1/2 


(24) 


I- 

The spatial coordinates of a domain wall section may x and x are unitary vector with the direction of the 

velocity and of the spatial directions of the domain wall. 
Assuming, without loss of generality that k = k&s, we 
have that 


be expressed, in this gauge, as 

X = Xo +6x -I- 6* 


uryx, (25) 

where xq is the (random) position of its center of mass, 

I 


000 = 


, . , sin ^fclx 3 ^^^/ 2 ^ sin ^fclx 3 ^^^/ 2 ^ 

dcryv^cos f k • X -I- vkrjxsj -- - , 


(26) 
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Qij — ©00 


v'^XiXj - (i - 


(1-^2) (x;«x:«+x;(^)x:(^)) 


(27) 


where the factor was included to compensate for 
the fact that we are only considering the real part of 
0 ^y(k, rj), and Xj = x- ej and Xj^'^ = x • ej are, respec¬ 
tively, the projections of vector x and x along the j-th 
spatial direction (defined by the unitary vector ej). 

For this choice of k, the Boltzmann integrator CMB- 
FAST — which is the basis of CMBACT - requires five 
components of the energy-momentum tensor [lOj. The 
first three are the scalar, vector and tensor components 
of the anisotropic stress, given by: 


where x is the position of the observer, n is the line 
of sight direction, and Tcmb is the average temperature 
of the CMB. The contributions of the different angular 
scales may be separated by doing a decomposition into 
spherical harmonics 

A(n) = ^ 

£m 

where are spherical harmonic functions. The angular 
power spectrum, is defined as 



2©s = 2©33 — ©11 — ©22 , 

(28) 


©V = ©13 , 

(29) 


©T = ©12 ■ 

(30) 

For single domain wall section, these yield 


205 = 

©00 ( 3 x 3 x 3 - 1 ) - 

(31) 
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(1 - v^) -f - 2 ) 

}. 

©y = 

©00 |w^XlX 3 - 

(32) 

- 

( 1 -u^) (x;«kY+x'/^)kY)}, 


0 J. = 

©00 |w^XiX 2 - 

(33) 


- (1 
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+ Xl X; 



The remainder components — the velocity field ©d 
( given by ©£) = ©03 for this choice of k), and the trace 
or isotropic pressure © = Qu — may be obtained from 
the covariant conservation of energy-momentum: 


©D = 


-2%d- 

a 

(©D- 0 Oo) -©00 


"O 

3 La 


2©s 


© = V ^©D — ©oo) — © 


00 ■ 


= 0, (34) 
(35) 


V. RESULTS AND DISCUSSION 

The CMB has a nearly perfect blackbody radiation 
spectrum [s^, with an approximately constant temper¬ 
ature accross the sky. For this reason, the basic CMB 
observables are its anisotropies, characterized by the tem¬ 
perature fluctuations 


A(x,n,ryo) 


|T(x, n,r]o) - TcmbI 


Tcmb 


(36) 


— 2£ I T/ , (38) 

m— — l 

where angled brackets represent an ensemble average. 

We calculate the CMB anisotropies generated by do¬ 
main wall networks by implementing the changes de¬ 
scribed in Sec. HVl to the publicly available CMBACT 
code (version 4.0). This code is based on CMBFAST 
which integrates the Einstein and Boltzmann equations 
using the line of sight method [Slj . CMBACT also com¬ 
putes the the cold dark matter (CDM) linear power spec¬ 
trum generated by active sources. 


P{k)=\S^{k)\ , (39) 

where d(k) is the Fourier transform of the density con¬ 
trast, 


w _ Pi^) - (P) 

^ 

p(x) is the matter density at a given position x, and {p) 
is its average value. 

Our results are found by averaging over 1000 dif¬ 
ferent realizations of a brownian domain wall network. 
We use the Planck 2015 cosmological parameters — 
= 0.0225 and = 0.1427 for the baryon 

and matter density parameters, respectively, and Hq — 
lOOh kms“^Mpc~^, with h = 0.6727, for the value of 
the Hubble parameter at the present time — and fix the 
domain wall tension to GaLo = 10“^ (where Lq is the 
characteristic lengthscale of the domain wall network at 
the present time). 

In Fig. m we plot the CDM linear power spectrum 
generated by domain wall networks along that of cosmic 
string networks for Gp = GaLo = 10“’^ (we use the un¬ 
altered CMBACT to obtain the power spectra for cosmic 
strings). As this figure illustrates, the matter power spec¬ 
trum generated by domain walls is strongly suppressed 
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FIG. 1: Comparison between the linear CDM power spectrum 
generated by domain wall (purple line) and cosmic string (red 
line) networks. We have averaged over 1000 realizations of 
both string and wall networks, and chose Gjj. = GaLo — 10“^. 



FIG. 2: Evolution of the linear CDM power spectrum gener¬ 
ated by domain wall networks. We include the CDM power 
spectra generated by domain walls until a = 10~^ (yellow 
line), a = 10~® (orange line), a = 10“^ (red line), a — 10“^ 
(blue line), and o = 1 (purple line). We have averaged 
over 1000 realizations of domain wall networks, and chose 
GaLo = 10-^ 


on small scales (or large fc), when compared to that of 
cosmic strings. Domain walls only become cosmologi- 
cally relevant at late cosmological times, when the cor¬ 
relation length of the network is large. For this reason, 
these networks should be expected to contribute mostly 
to the matter power spectrum on large scales. Fig. [5] 
— where the evolution of the power spectrum generated 
by domain wall networks is plotted — shows that this is 
indeed the case: the dominant contribution to the mat¬ 
ter power spectrum at low k is generated at late cosmo¬ 
logical times, while the large k contributions are mostly 
generated at earlier times (for a < 10“^). Note how¬ 
ever that, even on large scales, the CDM power spectrum 
generated by domain walls has an amplitude that is sig¬ 
nificantly smaller than that generated by cosmic strings. 
Domain walls have a correlation length that is slightly 
larger than that of cosmic strings at late times, which 
corresponds to a smaller energy density at the present 
time if G/i = GtrLo. This difference accounts for a factor 
of about 2 in the amplitude of the CDM power spec¬ 
trum. However, the observed difference is slightly larger 
because cosmic strings become cosmologically relevant 


earlier than domain walls, thus contributing significantly 
to the CDM power spectrum at large scales over a longer 
period of time. 

These differences in the shapes of the matter power 
spectra generated by domain walls and cosmic strings 
must necessarily translate into significant differences in 
the CMB signatures of these defects. In Fig. |31 we plot 
the TT, EE, TE and BB angular power spectra generated 
by domain wall networks along that generated by cosmic 
strings. First of all, the overall magnitude of the an¬ 
gular power spectra generated by domain wall networks 
on large angular scales (or small £) is lower than that 
of cosmic strings. Moreover, the strong suppression of 
the CDM power spectrum on small scales also results in 
a significantly reduced angular power spectra at large i 
when compared to that of cosmic string networks. In 
particular, the peak that is observed in the cosmic string 
temperature power spectrum at intermediate scales is ab¬ 
sent in the case domain walls. Unlike cosmic strings — 
which exhibit the aforementioned peak around i ^ 200 — 
domain walls contribute to the temperature power spec¬ 
trum mostly on large scales. For domain walls, then, the 
multipole modes around i = 2 have the highest constrain¬ 
ing power on their fractional contribution to the observed 
temperature power spectrum. Note that the uncertain¬ 
ties at low i are very large due to cosmic variance. For 
this reason, as discussed in [s^, observational data al¬ 
lows, on large angular scales, for a fractional contribu¬ 
tion of cosmic strings to the TT power spectrum around 
ten times larger than the upper limit obtained using the 
full data set (which is more constraining mainly due to 
the peak located at intermediate angular scales, where 
the observational uncertainties are very small). Current 
Planck data allows for a fractional contribution of cosmic 
strings of about 1-2% Em, and therefore one may 
conclude that the upper limit for the fractional contri¬ 
bution of cosmic strings in the TT channel using low 
I data only should be around 10-20%. One should ex¬ 
pect the maximum contribution of domain walls to the 
TT power spectrum allowed by observational data to be 
similar to that allowed for cosmic strings at large scales. 
By taking the conservative approach and assuming the 
fractional contribution of domain walls to temperature 
power spectrum at / = 2 to be around 20%, we obtain 
a constraint on the domain wall mass per unit area of 
a < 3.52 X 10“^ kgm“^, which corresponds to a con¬ 
straint on the domain wall-forming symmetry breaking 
scale of T] < 0.92 MeV. This estimate, despite being con¬ 
servative, is in good agreement with the Zel’dovich bound 
Q — which constrained r| to be smaller than 1 MeV — 
and even constitutes a slight improvement. Note how¬ 
ever that the large cosmic variance associated to the 
CMB anisotropies generated by standard domain wall 
networks on large cosmological scales must necessarily 
be dealt with in a more rigorous study. 

The aforementioned overall suppression of the temper¬ 
ature power spectra of domain walls (that becomes more 
accentuated with increasing £) is also observed in the 
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FIG. 3: Comparison between the angular power spectra generated by domain wall (purple lines) and by cosmic string (red 
lines) networks. From top to bottom, we plot the TT, EE, TE and BB power spectra, as a function of the multipole moment 
The left, middle and right panels represent the scalar, vector and tensor components, respectively. In the case of the TE 
power spectra, we chose to plot the absolute value and we use a dashed line to represent the parts in which Ct. is negative. We 
have averaged over 1000 realizations of both string and wall networks, and chose G/r = GaLo — 10~^. 


polarization power spectra, and it is slightly more ac¬ 
centuated on large scales in the EE and BB channels. 
However, since the upper bounds on GuLq are less strin¬ 
gent than those on G/i, the contribution of domain walls 
to the polarization power spectrum at large scales may in 
fact be significant. In Fig. |4l we plot the total BB power 
spectrum generated by domain wall networks and cosmic 
string networks that have the maximnm fractional con¬ 
tribution to the TT power spectrum allowed by current 
observational data. We chose cr = 3.52x10“® kgm“^ (or 
equivalently GctLq = 5.6 x 10“®) for the domain wall net¬ 
works, and the weakest constraint on cosmic string ten¬ 


sion obtained using Planck data [s^, G/r = 2.4 x 10“^ 
(which corresponds to the constraint on Abelian-Higgs 
strings). This figure shows that the B-mode polarization 
signal generated by domain wall networks that have a 
fractional contribution to the temperature power spec¬ 
trum that is close to the allowed observational limit may 
dominate over that generated by cosmic strings (and the 
same is also true for the EE signal). The TT, EE and 
TE signals generated by topological defects are expected 
to be subdominant, when compared to that generated 
by primordial fluctuations. However, their contribution 
in the B-mode channel may be dominant, due to the 

























FIG. 4: Comparison between the total BB power spectra gen¬ 
erated by domain wall (purple line) and cosmic string (red 
line) networks with the maximum tension allowed by current 
observational data. For domain walls, we have used the up¬ 
per limit obtained in this paper GctLq = 5.6 x 10“®, while for 
cosmic string we took the upper limit obtained using Planck 
data Gfi = 2.4 x 10“^ [3^ . We have averaged over 1000 real¬ 
izations of both string and wall networks. 

presence of a vector component contribution that is ab¬ 
sent in inflationary scenarios. Our results seem to indi¬ 
cate, therefore, that domain walls could be the dominant 
contributor in the BB channel for low multipole modes. 
Moreover, the fact that domain walls and cosmic strings 
have spectra with different shapes and that these defects 
contribute mostly at different scales should, in principle, 
allow to distinguish between these contributions if a sig¬ 
nal is detected. 

As is the case for the matter power spectrum, the dom¬ 
inant contributions to the CMB anisotropies are gener¬ 
ated in the matter era. In Fig. [SJ we plot the evolu¬ 
tion of the TT, EE, TE and BB power spectra generated 
by a domain wall network. These plots show that do¬ 
main walls contribute to the temperature anisotropies on 
progressively larger angular scales (or lower multipoles) 
as time progresses and their correlation length increases 
(as was the case for the CDM power spectrum). For 
the polarization power spectra, however, the picture that 
emerges is slightly different. Polarization may be created 
in two narrow windows in the history of the universe: 
very close to recombination and during the reionization 
epoch (because it requires the presence of a temperature 
quadrupole and the universe to be ionized). Fig. [S] illus¬ 
trates this fact: most of the contributions to the small 
scale peak are generated around the last scattering epoch 
(while a ^ 10“^), while the large scale (dominant) peak 
is mostly generated at more recent cosmological times, 
around the time the scale factor was a ^ 10“^. 


VI. CONCLUSIONS 

In this paper, we have expanded CMBACT to allow 
for the numerical computation of the CMB anisotropies 


generated by domain wall networks. This was done by 
adapting the USM for cosmic strings in order to also 
allow for the description domain wall networks, and by 
implementing the necessary changes on CMBACT (ver¬ 
sion 4.0). Note that, within this framework, the VOS 
model is used to set the correlation length and average 
velocity of the domain wall sections. This phenomenolog¬ 
ical model includes the essential aspects of domain wall 
dynamics, and it has the advantage of having two free 
parameters that may be used to calibrate the model in 
order to describe different types of domain wall networks. 
Our approach, thus, has the advantage of not being lim¬ 
ited to a specific domain wall scenario. 

We have also used this tool to study the angular power 
spectra generated by standard domain wall networks, in 
both the temperature and polarization channels, and to 
derive a conservative constraint on the energy scale of 
formation of the domain wall network of r] < 0.92 MeV. 
Note, however, that this only applies to standard domain 
wall networks. Other domain wall scenarios — such as 
networks whose dynamics is friction-dominated or do¬ 
main wall networks with junctions — would be charac¬ 
terized by a smaller correlation length. L determines 
in which multipole mode the networks contribute domi¬ 
nantly to the temperature power spectrum. One should 
therefore expect constraints on non-standard networks to 
be stronger than the one obtained here, since they would 
be constrained by higher multipole modes which have 
smaller observational uncertainties and which are less af¬ 
fected by cosmic variance. The study of these different 
scenarios will be the subject of future work. 

Note: While this manuscript was in preparation, re¬ 
sults related to the ones presented here appeared in Ref. 
[S^. The shape of the power spectra obtained in this 
work seems to be in good qualitative agreement with our 
results and the constraint on t| presented by the authors 
is very similar to the one we have obtained. However, 
their constraint on a is four orders of magnitude smaller 
than the constraint we have obtained and it seems to be 
inconsistent with their upper limit of rj. 
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SCALAR VECTOR TENSOR 




FIG. 5: Evolution of the angular power spectrum generated by domain wall networks. From top to bottom, we plot the TT, 
EE, TE and BB power spectra, as a function of the multipole moment £. The left, middle and right panels represent the scalar, 
vector and tensor components, respectively. In each of the plots we include the angular power spectra generated by domain 
walls until a — 10“'* (yellow line), a = 10“® (orange line), a = 10“^ (red line), a = 10“^ (blue line), and a — \ (purple line). In 
the case of the TE power spectra, we chose to plot the absolute value and we use a dashed line to represent the parts in which 
Cl is negative. We have averaged over 1000 realizations domain wall networks, and chose GgLq = 10“^. 
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